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Abstract- Conventional  screening  for  Fetal  Alcohol  Syndrome 
in  children  involves  time-consuming  and  intrusive  facial 
measurements.  We  are  developing  an  alternative  method  that 
will  allow  automatic  measurement  of  the  relevant  points  from 
a  pair  of  stereo  photographs.  This  paper  describes  a  method  for 
automatic  eye  extraction  from  such  photographs,  using  genetic 
algorithms  to  match  eye  templates  to  face  images. 
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I.  INTRODUCTION 

Fetal  Alcohol  Syndrome  (FAS)  is  a  consequence  of  the 
teratogenic  effects  of  excessive  maternal  alcohol  ingestion  on 
the  fetus  in  utero ,  and  is  the  most  common  preventable  cause 
of  mental  retardation  worldwide  [1],  [2].  FAS  is  characterized 
by  growth  retardation,  central  nervous  system  abnormalities, 
facial  dysmorphology,  and  other  malformations.  Atrophy  and 
flattening  of  the  frontal  lobes  of  the  brain  are  thought  to  cause 
the  characteristic  facial  features,  mainly  contraction  of  the 
middle  third  of  the  face. 

There  are  no  biological  markers  for  FAS.  Diagnosis  requires 
a  positive  maternal  history  of  heavy  alcohol  ingestion  during 
pregnancy  in  addition  to  the  presence  of  the  characteristic 
features,  confirmed  by  specialists  with  experience  and 
expertise  in  dysmorphology.  Conventional  assessment  of  facial 
dysmorphology  includes  time-consuming  and  intrusive 
distance  measurements  on  the  face.  The  relevant  points  and 
distances  are  depicted  in  Fig.  1.  We  are  developing  an 
alternative  method  that  will  allow  automatic  measurement  of 
the  relevant  points  from  a  pair  of  stereo  photographs  of  6  to  7- 
year  old  children.  This  paper  describes  the  use  of  genetic 
algorithms  (GAs)  to  find  eye  features  in  these  photographs. 
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Figure  1 .  Distance  measurements  for  FAS  screening 


A.  Image  Acquisition 

Two  photographs  are  taken  simultaneously  with  digital 
cameras  mounted  about  40  cm  apart  at  a  fixed  distance  from 
a  reference  frame  within  which  the  children  place  their 
heads  (Fig.  2).  Markers  on  the  reference  frame  allow  the 
images  to  be  calibrated  and  placed  in  a  three-dimensional 
co-ordinate  frame.  If  the  eyes  are  located  and  the  relevant 
points  identified  on  each  of  the  two  images,  the  required 
distances  can  be  calculated. 


Figure  2.  Stereo  photographs  using  reference  frame 
B.  Genetic  Algorithms 

GAs  are  adaptive  search  procedures  that  perform  function 
optimization  based  on  the  principles  of  evolution  by  natural 
selection  [6].  They  are  capable  of  locating  good 
approximate  solutions  in  large  search  spaces  without 
becoming  trapped  in  local  minima. 

The  function  being  optimized  by  the  GA  is  known  as  the 
objective  function  or  fitness  function,  and  its  range  is  called 
the  search  space.  The  parameters  of  the  optimization 
problem  are  represented  by  binary  strings  called  genes  that 
are  concatenated  into  chromosomes  representing  potential 
solutions.  If  Yi  denotes  a  gene  representing  parameter  i,  then 
the  chromosome  c  =  Yi>Y2>--->Yn>  and  a  population  P  =  (ch 
c2, ...,  ck},  where  k  denotes  the  population  size. 

The  initial  population  of  chromosomes  with  which  the 
search  begins,  is  generated  randomly.  Each  member  of  the 
population  is  evaluated  with  the  fitness  function  and 
assigned  a  fitness  level.  The  population  is  then  iteratively 
updated  through  a  reproduction  strategy  that  favours  the 
fittest  chromosomes,  until  the  fitness  function  is  maximized 
or  the  desired  number  of  generations  has  been  reached. 

GAs  are  employed  here  to  detect  eyes  in  an  image  by 
finding  the  optimal  match  firstly  between  a  circular  iris 
template  and  the  image,  and  then  between  the  sclera 
template,  described  by  cubic  splines,  and  the  image. 
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C.  Eye  Templates 

Yuille  et  al.  [3]  proposed  and  Xie  et  al.  [4]  modified  a 
method  of  eye  feature  extraction  based  on  deformable 
templates  in  which  the  iris  is  modelled  by  a  circle  and  the 
eyelids  by  2  parabolic  sections.  Esme  et  al.  [5]  used  genetic 
algorithms  as  a  search  strategy  to  find  eye  contours  modelled 
by  parabolic  templates. 

We  propose  an  algorithm  in  which  the  iris  is  modelled  by  a 
circle  as  suggested  in  [3]  and  [4],  but  the  eye  by  a  set  of  cubic 
splines  (one  representing  the  upper  and  one  the  lower  eyelid). 
Each  cubic  spline  is  described  by  4  control  points,  and  the  end 
points  are  common.  The  eye  template  is  shown  in  Fig.  3. 


Figure  3.  Eye  template.  The  circle  is  specified  by  its  radius  and  the  co¬ 
ordinates  of  its  centre.  The  arrows  indicate  cubic  spline  control  points. 


normalized  to  the  range  0  to  255.  I  denotes  the  original 
image,  and  u  denotes  the  unit  step  function. 

The  peak  map  (Ip)  selects  bright  regions  in  the  image, 
such  as  the  sclera,  in  a  similar  manner:  it  assigns  the  value 
‘O’  to  those  pixels  less  than  the  window  average,  and 
assigns  to  those  pixels  greater  than  the  window  average  a 
value  proportional  to  the  difference  between  the  pixel  and 
the  average. 

E.  Matching  the  Iris 

1)  Parameters'.  The  parameters  being  optimized  by  the 
GA  and  which  are  therefore  encoded  in  the  chromosomes 
are  the  radius  and  centre  (x-  and  y-co-ordinates)  of  the 
circle. 

2)  Initial  population :  An  initial  random  population  of 
chromosomes  is  generated. 

3)  Fitness  function:  The  iris  appears  as  a  bright  circular 
area  in  the  valley  image  Iv.  The  circular  template  is  filled  in 
so  that  it  becomes  a  disc  Id.  The  correlation  coefficient  r  of 
this  disc  with  the  valley  image  is  the  fitness  function  f  as 
shown  in  (3)  and  (4).  A  and  B  in  (3)  denote  any  two 

images,  while  A  and  B  denote  the  mean  pixel  intensity  in 
these  images. 


D.  Pre-processing 

Since  the  face  is  always  photographed  in  a  fixed  position,  an 
image  window  in  which  each  eye  can  be  expected  to  appear 
can  be  specified  in  order  to  reduce  the  search  space.  Each  eye 
is  treated  separately. 

Yuille  et  al.  [3]  and  Xie  et  al.  [4]  used  search  strategies 
based  on  steepest  descent  and  the  Levenberg-Marquardt 
method  respectively,  operating  on  complex  energy  functions 
consisting  of  terms  due  to  valley,  edge,  peak,  image  and 
internal  potentials,  to  find  the  desired  points  in  the  eye.  Esme 
et  al.  [5]  used  a  fitness  function  based  on  fused  peak,  valley 
and  edge  maps.  Our  fitness  function  is  based  only  on  peak  and 
valley  maps. 

The  peak  and  valley  maps  [5]  of  the  image  are  used  to  fit  the 
templates  to  the  image.  They  are  calculated  as  shown  in  (1) 
and  (2). 


T(  X  u(lw(x,y)-l(x,y))xlw(x,y)-l(x,y) 
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(2) 


The  valley  map  (/v)  selects  dark  regions  in  the  image,  such 
as  the  iris.  It  compares  pixels  in  an  image  window  (Iw)  with  the 
window  average,  assigns  the  value  ‘0’  to  those  pixels  greater 
than  the  window  average,  and  assigns  to  those  pixels  less  than 
the  window  average  a  value  proportional  to  the  difference 
between  the  pixel  and  the  average.  The  new  pixel  values  are 
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(4) 


4)  Reproduction:  One  third  of  the  chromosomes  in  each 
new  generation  are  those  with  the  highest  fitness  values 
from  the  previous  generation.  Another  third  of  the  new 
generation  is  obtained  by  performing  simple  crossover 
between  randomly  selected  pairs  of  strings  from  the 
previous  generation  with  no  weighting  based  on  string 
fitness,  and  the  final  third  consists  of  new  randomly 
generated  chromosomes.  The  population  size  remains 
constant. 

F.  Matching  the  Sclera 

1 )  Parameters:  Once  the  iris  has  been  localized,  a  smaller 
window,  within  which  the  sclera  is  expected  to  be  located, 
is  defined  around  the  iris  centre.  The  parameters  being 
optimized  are  the  locations  (x-  and  y-co-ordinates)  of  the  6 
cubic  spline  control  points  that  specify  the  eye  template. 

2)  Initial  Population:  The  initial  population  is  generated 
randomly  as  with  the  iris. 

3)  Fitness  function:  The  sclera  appears  as  a  bright  region 
in  the  peak  image  Ip.  The  area  within  the  2  cubic  splines  that 
form  the  eye  template  is  filled  in  and  the  iris  is  removed 
from  it  so  that  only  the  sclera  remains.  This  template  is 
multiplied  point- wise  with  a  gaussian  of  the  same  size  so  as 


to  simulate  the  effect  found  in  images  of  the  eye  and  the  peak 
map  where  the  intensity  of  the  sclera  gradually  decreases 
towards  the  eye  corners. 

The  valley  image,  in  which  the  sclera  appears  dark,  is  also 
used  in  the  fitness  function,  which  is  shown  in  (5).  Is  denotes 
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Figure  8.  Results  for  four  different  sets  of  eyes. 


IV.  CONCLUSION 

We  have  presented  an  algorithm  for  eye  extraction  that 
we  intend  using  as  part  of  a  method  of  screening  children 
for  the  facial  features  characteristic  of  fetal  alcohol 
syndrome.  The  accuracy  at  which  the  desired  points  on  the 
eye  (corners,  iris  centre)  can  be  located  is  yet  to  be 
determined.  Future  work  includes  comparison  of  the 
performance  of  this  algorithm  with  manual  identification  of 
the  relevant  points. 

ACKNOWLEDGMENT 

This  work  was  conducted  with  the  aid  of  a  grant  from  the 
South  African  Medical  Research  Council. 

REFERENCES 

[1]  D.Viljoen,  “Fetal  alcohol  syndrome”,  South  African  Medical  Journal , 
vol.  89(9),  pp.  958-960,  1999. 

[2]  J.  Croxford,  D.  Viljoen,  “Alcohol  consumption  by  pregnant  women  in 
the  Western  Cape,”  South  African  Medical  Journal,  vol.  89(9),  pp.  962- 
965,  1999. 

[3]  A.L.  Yuille,  P.W.  Hallinan,  D.S.  Cohen,  “Feature  extraction  from  faces 
using  deformable  templates,”  Int.  J.  Computer  Vision,  vol.  8(2)  pp.  99-111, 
1992. 

[4]  X.  Xie,  R.  Sudharkar,  H.  Zhuang,  “On  improving  eye  feature 
extraction  using  deformable  templates,”  Pattern  Recognition,  vol.  27(6), 
pp.  791-799,  1994. 

[5]  B.  Esme,  B.  Sanur,  E.  Anarum,  “Facial  feature  extraction  using  genetic 
algorithms,”  Signal  Processing  VIII,  Theories  and  Applications, 
Proceedings  of  EUSIPCO-96,  1996. 

[6]  D.  E.  Goldberg,  Genetic  Algorithms  in  Search,  Optimization,  and 
Machine  Learning,  New  York:Addison  Wesley,  1989. 


